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INTEGRAL EQUATIONS, FRACTIONAL CALCULUS AND 
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Abstract. We present an extension of a previously developed method em- 
ploying the formalism of the fractional derivatives to solve new classes of in- 
tegral equations. This method uses different forms of integral operators that 
generalizes the exponential shift operator. 



1. Introduction 



(-H In a previous note [T we have obtained a solution of the Lamb integral equation 

pi I [21 [3] using a formalism based on the fractional derivatives. Our result coincides 

with the solution provided, without any proof, by Bateman jS^. A generalization of 

our method has been proposed in [4, , where fractional derivatives are used to solve 

a more complicated version of the Lamb equation, 
y In this note we discuss a further extension of the technique put forward in ^ 

and we will develop a more general procedure to treat integral equations whose 
I solution requires fractional forms of differential operators other than the ordinary 

^ derivative. To this aim we remind that the exponential operator 

^ (1) EiX) ^ exp{Xxd^} , 

ly-v is a dilatation operator whose action on a given function g(x) is (see ref. [S]) 

^ (2) E{X)gix) = g{e^x). 

^D Let us now consider the operator {x dx)~^ with 3fi i^ > 0. The use of the following 

^"^ property of the Laplace transform 

' -• —sa ^v — l 



C^ allows to write 

(4) (^a^)-^_y^ dsE{^s)s''-\ 

and, according to the previous discussion, we find [6] 

(5) (^a^)-g(^) = _y^ AsE{-s)g{x)s''-' 



which is essentially the generalization of the Riemann-Liouville integral represen- 
tation for non-integer negative powers of the operator x S^; [7] . 
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2. Integral equations and fractional differential operators 

The procedure outlined in sec. [T] can be generalized to other differential forms. 
The first example we consider is the following modified form of the Lamb-Bateman 
equation 

/■oo 

(6) / dyu{e~y\) = fix), 

Jo 

where u{x) is the function to be determined, and f{x) is a continuous, free from sin- 
gularities, known function. According to the properties of the dilatation operator, 
we can cast eq. ^ in the form 

/>oo 

(7) / dye-y"^''-u{x)^ fix), 

Jo 

and, treating the derivative operator as a generic constant, the evaluation of the 
the gaussian integral yields 
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(8) uix) = —- ^xdxfix). 

By rewriting the operator \Jxdx as ixd^) ix dx)~^''^ , and taking into account eq. 
@, we gelQ 
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whose correctness as solution of eq. ^ has been checked numerically. 

As a further example of integral equation which does not involve fractional op- 
erators but requires the use of the dilatation operator, we consider the following 
problem 

(10) / dye-yuiyx) = fix), 

Jo 

where uix) is the unknown function. In this case we get 

/•oo 

(11) f{x)= dye-yy^^'-uix) ^Vixd^ + l)uix), 



and, thus 

(12) uix) = {Vixd, + l)}-^ fix). 

In the case fix) = e~^ , we obtain the solution in terms of 0-th order Bessel 
functions of first kind, nameljrl 

(13) uix) ^ V LJ- {Tixd, + l)}-' x" 



oo 



( — \^ 



Let us remark that, as a consequence of the relation [/(O), O] = 0, with / generic function 
that admits a power series expansion, one has: (xdx) {xdx)~^'^ = {xdx)~^'^ (xdx)- 

Note that from {xdx)x" = nx" we obtain the following identity g{xdx)x" = g{n)x". 
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In general, if /(x) is any function which can be expanded as 

oo 

(14) /(x) = ^^x", 

the solution writes 



n=0 



(15) ,(,) = ^^,V 

It is interesting to note that in the case in which the equation is of the type 

/"OO 

(16) / Ay^ryuiiTx) = e"^ , 

the relevant solution can be written in terms of Bessel-like functions as follows 

(17) u(x) = y" M- {r(TOxa:, + i)}"^ X" 

^ — ' n\ 
n=a 

°° (' — "l" 

= E u V ^" = w^o(xH. 

^-^ n\ (mn]\ 

n=0 ^ ' 



where 



oo ( \k 

(18) W^{x\m) = Y. -' - ^' 



fc=0 



fc!(m fc + n)\ 



is the Bessel- Wright function of order n [S] . The same result holds for the integer 
m replaced by any real /i . 

We go back to fractional operational calculus by considering the following exam- 
ple of integral equation 

/■oo 

(19) / d2/u(v/x2 + 2y2) = f{x) 

Jo 

which is apparently more complicated than the original Lamb-Bateman equation, 
but, as will be shown later, in spite of its different form, it represents the same 
mathematical problem. Also in this case the solution can be obtained using the 
previously outlined shift operator technique. We first recall the identity [6] 

(20) exp |a ^ a,| g{x) = g{Vx^ + 2X) 

and therefore, as a straightforward consequence of our procedure, we find 



(21) u{x) = -^J--d^f{x). 

^Ji: V X 

The use of the same logical steps leading to eq. (|9| yields the following expression 
for the solution of eq. ( 19 ) 

(22) .,.,.-^rd,-M(|) ^ 
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3. Concluding remarks 

The examples discussed so far show that the method we have proposed is fairly 
powerful and is amenable of useful generalizations. A further example supporting 
the usefulness of our technique is offered by the calculation of the integral 

(23) F{x-v)=j^ disin^^^^^J (^v > 

The shift operator method outlined in the introductory remarks allows to write this 
integral as follows 

POO / 1 \ V Xd:c 

(24) F{x;u) = Oy sinx with 6,. = / dt ( ^ 

J a V 1 + ^ 

The operator 0^ can be defined explicitly by using the following result [5] 

1 _ V^ r(/i-i/2) 
2 v{^,) ' 



(25) 



'^^'{i + y^Y 



that specified to the case /x = vxdx, yields 

y^T{iyxdx-l/2) 
2 



(26) 



F{x;,y) 



sma; 



r{vxdx) 

By substituting to sin a; its Taylor expansion, and assuming that the sum and the 
operator O^ can be interchanged, we get 



y OO 

(27) F{x;,^) = ^r{:,xdx-l/2){ri,.xdx)}-'J2 



(~r 



^E 



(-)" r((2n+l)i^-l/2) 



(2n + l)! 



,2n+l 



„2n+l 



2 ^^ (2n+l)! r((2n + l)^) 



We have checked the correctness of this result performing a direct numerical inte- 
gration of eq. ( p3| . The behavior of the functions F{x; v) for different values of v 
is shown in Fig. [T| 

The procedure employed for the evaluation of the integral ( 23 1 is quite general 
and can be extended to any integral of the type 

/"OO 

(28) G{x) = / Atf[xg[t)) 

Jo 

In fact, writing this integral as follows 



(29) 

we can calculate 

(30) 

and, thus 
(31) 



G{x) 



dt [.g(t)]"'^ fix) 



Oif,) = / dt [g{t)r , 

^0 



Gix) = Oixdx)fix). 

The use of the shift operator method allows also the evaluation of the following 
integral 



(32) 



/■CSO 

I{x) = / dtf{^x^-2git)), 
Jo 
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Figure 1. The function F{x:v) for v — 1.5 (dot) and v — 4.1 (solid). 



which, on account of cq. ([20| , takes the following operational form 



(33) 
where 

(34) 



1 



I{x) ^Q\-dA fix) 



Jo 



A complete theory of the generalized shift operators has been presented in ref. 
[5], where, among the other things, the following identity has been proved 
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(35) exp {A q{x) d,} g{x) = g [p-^ (A + F{x))\ F{x) 



de 



m 



which is a generalization of the ordinary shift operator (the examples discussed 
before are just a particular case of this formula). According to this identity and 
to the methods here discussed, it is easy to show that the solution of the following 
integral equation 



(36) 

can be written as 

(37) 



dyu[F-^{-y^ + F{x))\ = f{x) 



u{x) 



\/q(x)Kf{x) 



In a forthcoming investigation we will discuss the conditions under which the pre- 
vious solution holds and some examples of applications of the obtained identities. 
Before closing this paper we consider the following integral equation 



(38) 



dyu 
\l + xy 



fix) 



By using eq. (35), it is easy to show thai 
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(39) " ( ^i— ) ^ ex.p {-y x"^ dx} u{x) 

\l + xy J 

and, therefore, for the solution of cq. ( [38| we obtain 

(40) u{x) = "^'^" ,„. fix) 

In this paper we have studied the properties of differential operators like x dx 
and their fractional generalizations. It is well known that, for n integer different 
from zero, the following identity holds [9] 

n 

(41) (xa,)" = ^5(n,A:)x'=9,\ 

A:=0 

where 

(42) S{n,k)^l-^Y^{-f-^(^\r 

are the Stirling number of second kind. The extension of this equation to the non- 
integer powers of x d^ does not exist. We make the conjecture that in this case eq. 



(41 1 can be modified as follows 



oo 



(43) [xdxY - y" S{k,n)x^dl (0 < i/ < 1) 



fc=0 



The numerical check confirms the validity of our conjecture, but we did not succeed 



in getting a rigorous proof for it. Eq. (43), as well as other identities presented in 



this paper, will be more throughly discussed in a more general forthcoming paper. 



This relation can also be obtained setting x = 1/^ (with xy ^ —1). 
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